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ABS 'TRACT 

The Transonic Turbine Test Rig at the Naval 
Postgraduate School, Monterey, California, uses a super- 
sonic jet exhauster to lower the turbine exit static 
pressure and thus increase pressure ratios available for 
turbine testing. Exhauster operation results in inter- 
action between two co-flowing streams of air. The primary 
air is supersonic and the secondary air, subsonic. This 
thesis is an experimental study of the mixing process that 
occurs between the two streams in a cylindrical mixing 
chamber. In addition, a theoretical model of the exhauster 
performance based on polytropic efficiencies and conserva- 
tion of momentum was programmed, and the predicted values 
were compared to test results. 
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1. Introduction 



Exhausters are being incorporated into the design of 
many systems because of the simplicity of their design and 
the absence of moving parts. Air augmentation of rocket 
thrust, torpedo drag reduction innovations, use as a jet 
pump in aircraft for fuel transfer or boundary layer con- 
trol are significant applications. 

The Transonic Turbine Test Rig (TTTR) installed at the 
Propulsion Laboratory of the Naval Postgraduate School uses 
a supersonic exhauster system to increase the range of tur- 
bine pressxire ratios available for testing. The test tur- 
bine is enclosed in a hood. A compressor can provide a 
maximxim turbine inlet total pressure of 43 psia giving a 

maximiam turbine pressure ratio (Po. /P_ ) of about 3 

^inlet exit 

when exhausting to the atmosphere. Since the turbine uses 
only part of the mass flow furnished by the compressor, the 
remainder can be used to drive the exhauster which produces 
a vacuum in the hood and thereby increases the turbine 
pressure ratio. 

Very little information is available on the mixing 
phenomena encountered with supersonic exhausters. 
Abramovich^ provides extensive information for the case of 

an inconpressible free jet and associated turbulent mixingl 
He includes test work accomplished, and theory developed, 

^Abramovich, G. N. , The Theory of Turbulent Jets 
(Cambridge, MassachusettsT The kiT iPress, Massachusetts 
Institute of Technology, 1963) , Part I, "Turbulent Jets 
in Incompressible Fluid.” 
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by Prandtl and Tollmien. The effect of compressibility on 
txirbulence is very complicated however, and simplifica- 
tions of the problem must be introduced to develop approx- 
imate solutions. 

This thesis is a study of the supersonic exhauster 
system used with the TTTR. Objectives of the investigation 
were to determine the hood total pressure variation with 
drive nozzle total pressure for a given secondary flow 
rate, and to study the mixing and turbulent interchange that 
occurs between the two co-flowing streams in the mixing 
pipe. Main topics include: traverse probe design, experi- 

mental results of the ejector test series, a computer pre- 
diction program of exhauster performance based on poly- 
tropic processes and conservation of momentum, empirical 
determination of polytropic efficiencies, and an attempt to 
determine from test data the losses due to mixing and tur- 
bulent interchange between the two streams of air. 

A summary of the general properties of turbulent jets, 
including terminology used to describe the different mix- 
ing regions that occxir with restricted jet operation and 
co-flowing streams, is given in Appendix D. The terms 
ejector and exhauster are used interchangeably throughout 
this thesis. 

The author is most appreciative of the guidance, coun- 
sel and encouragement provided by Professor M. H. Vavra, 
the laboratory assistance given by Mr. J. E. Hammer, and 
the probe fabrication work of Mr. T. B. Dunton. 
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2. Installation. 



The TTTR installation is shown schematically in Fig. 1, 
and Fig. 2 is a photograph of the turbine hood and exhauster 
assembly. High pressure air is supplied by an Allis- 
Chalmers VA 312 compressor and enters Tank 1. Air to the 
turbine passes from Tank 1 through a calibrated flow mea- 
surement nozzle and enters Tank 2. The amount of air to the 
turbine is controlled by a manual valve; an electrically 
operated exhauster inlet valve controls the drive nozzle 
total pressure (PTN) . More detailed descriptions of the 
TTTR and its operating performance are given by Eckert^ and 

3 

by Commons. 

The exhauster section of the TTTR is shown in greater 
detail in Fig. 3. Important features are the arrangement 
of twelve mixing tube static pressure taps (M2 - Ml3) , the 
Kiel probe to measure drive nozzle total temperature (TTN) 
and total pressure, and the secondary air static pressure 
taps just prior to entry into the mixing pipe. Secondary 
air is defined as the air that enters the mixing pipe from 
the turbine hood. To obtain the mixing pipe static pres- 
sures, four taps are equally spaced circumferentially at 



^Eckert, R. H. , Performance Analysis and Initial Tests 
of a Transonic Turbine Test Rig (Aeronautical Engineer's 
Thesis, U. S. Naval Postgraduate School, Monterey, 
California, May 1966), Section 2, "Installation." 

^Commons, P. M. , Instrumentation of the Transonic 
Turbine Test Rig to Determine the Performance of Turbine 
Inlet Guide Vanes Through the Application of the Momentum 
and Moment of Momentum Equations (Master's Thesis, U. S. 
Naval Postgraduate School, Monterey, California, 1967) , 
Section 2, "Installation." 
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each location and cominoned to obtain an average reading. 
Secondary air static pressxires are measured by three taps 
at each location, equally spaced circumferentially and 
commoned. 

One section of the mixing pipe was modified by cutting 
a 3 inch by 1/4 inch slit at 3 locations to provide ingress 
and egress of a traverse probe. The slits were designed to 
place the traverse probe total pressure opening at the same 
distance from the drive nozzle exit as the mixing pipe wall 
static pressure; that is, in the plane of the wall static 
pressure taps for a given axial location. Seals were also 
fabricated to prevent leakages, with or without the tra- 
verse probe in position. 

Secondary air total temperature (TTD) is measured with 
a modified Kiel probe. Hood total pressure (PTD) is mea- 
sured by eight static pressure ports that are commoned to 
obtain an average value. The secondary air is essentially 
that which passes from Tank 2 to the turbine. It is aug- 
mented by a small amount that enters the shaft labyrinth 
seal of the hood. Dimensions of the shaft labyrinth shown 

4 

in Fig. 4 differ from those given by Eckert due to increas- 
ing shaft clearance from 0.005 inch to 0.020 inch. A 
thermocouple is placed at the labyrinth entrance to measure 
the total temperature (TTC) of the leakage air. 

Drive nozzle I, the exhauster drive nozzle designed for 
the TTTR, is shown in Fig. 5 with detailed measurements. 

4 

Eckert, op. cit. , Fig. 4, p. 25. 
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Drive nozzle II, a smaller nozzle with the same area ratio 
(Athroat^^exit^ , was designed and built to investigate the 
effects of different primary and secondary flow areas, and 
smaller primary flow rate, on the mixing process and on hood 
total pressure. Details of Drive nozzle II are shown in 
Fig. 6. Design calculations were made to insure a uniformly 
decreasing secondary flow area for the new nozzle. 

All temperatures are measured using standard Iron - 
Constantan or Cromel - Alumel thermocouples with an ice bath 
as reference. Differential pressures across the flow mea- 
siaring nozzle are measured in inches of water on a mano- 
meter board with 0.1 inch gradations. Mixing tube static 
pressures are measured against atmospheric pressure with a 
Texas Instruments Fused Quartz Pressure Gage. All other 
pressures are measured in inches of mercury against atmos- 
pheric pressure as reference using 96 inch long manometer 
tubes with 0.1 inch gradations. A check of static pressxire 
readings between the Texas Instrument gage and the mercury 
manometers produced values within 0.5 per cent of each 
other. 

3. Probe Design. 

In order to study the mixing of two co-flowing streams, 

velocity traverses were desired at a number of axial loca- 

5 

tions. Naviaux attempted to obtain this information but 

^Naviaux, J. C., Transonic Turbine Test Rig Exhauster 
^stem Tests and Tests of a Reaction Turbine (Master's 
•nxesis, U. S. Naval Postgraduate School, Monterey, 
California, December 196b), Section 4, "Test Results." 
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experienced difficulty in obtaining reliable data. During 
his tests the probe had excessive vibration, the static 
pressure probe used had the pressure ports only 2 diameters 
away from the probe tip, and a thermocouple was not incor- 
porated, hence, the differences in stagnation temperatures 
between the two streams could not be measured. 

For the ejector tests, the four probes of Fig. 7 were 
designed and built to obtain traverse information; two were 
pressure probes with a thermocouple and two were hot wire 
probes. Probe I was the only successful instrument and all 
traverse information was obtained with it. 

Probe I is shown with its significant dimensions in 

g 

Pig. 8. Criteria presented by Gettelman and Krause con- 
cerning static port location for maximum accuracy were 
strictly adhered to. Best results are obtained when the 
static ports are 4-6 sensing element diameters away from 
the tip, and greater than 8 support diameters away from the 
support. (These minimum criteria produce static pressure 
readings within 5.0 per cent of the velocity head at 
Mach 0.9 and 3.5 per cent at Mach 0.5. Mr. E. Keener of 
NASA Ames Laboratories assured the author that the error 
decreases in supersonic conditions. 



^Gettelman, C. C. and Krause, L. N., Considerations 
Entering into the Selection of Probes for Pressure 
Measurement in Jet Engines (ISA Proceedings Vol. 7, 
Aeronautic Instrumentation - Session II), pp. 134-137. 
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The blockage of probe I is 4.3 per cent of the mixing 
7 

pipe area. Tests show that a blockage effect based on ratio 
of support diameter to duct or jet diameter can induce an 
error as large as 8.0 per cent for conditions and dimensions 
used in the ejector tests. The error, however, causes a 
measurement that is too large in a circular duct and a 

measurement that is too small in a jet. This difference in 
direction of the measurement error probably leads to a 
partial cancellation of the blockage error at mixing stations 
near the drive nozzle for the exhauster system. The exact 
effect and variations with axial position were not evaluated. 

The thermocouple of probe I was placed in the support 
0.5 inch above the pressure element. This, position was 
chosen to avoid a second extended sensing arm which could 
induce errors due to probe interactions » The thermocouple 
is exposed to the stream by a 0.160 inch diameter hole. A 
0.040 inch diameter hole in the after side provides air 
circulation. Total temperature measurements at mixing pipe 
station Ml were less than 2.0 per cent different than those 
obtained for the primary and secondary flow before entry to 
the mixing pipe. 

The probe vibration problem was eliminated by instal- 
lation of two 4 inch x 6 inch box supports contoured to fit 
the outside surface of the mixing pipe. Each box was 

7 

Krause, L. N. and Gettelman, C. Co, Effect of Inter- 
action Among Probes, Supports, Ducts, Walls and Jet 
Boundaries on Pressure Measurements in Ducts and Jets (ISA 
Proceedings Vol. 7, Aeronautic Instrumentation - Session II), 
Fig. 4, p. 140. 
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fitted with a 3 inch teflon bushing to provide probe support 
over a long distance. This eliminated the pivot point about 
which a probe would tend to oscillate. The assembly is 
shown in detail in Pig. 9 and pictured in Fig. 10. The 
design proved very successful. 

The mixing process is primarily a fvmction of axial 
location in the mixing pipe. Probe I has a 3 inch axial 
distance between the position where the total pressure is 
taken and the thermocouple. Because significant gradients 
in total temperature exist, particularly near the drive 
nozzle, probe II was designed to obtain the temperature 
nearer to the total pressure opening. An exposed thermo- 
couple was placed in a position 1/2 inch aft of the static 
pressure ports. 

Testing of probe II produced e^qjected results in the 
stibsonic secondary flow. In the supersonic region of the 
primary flow, the test data were irregular both for static 
pressures and temperatures. However, lower temperatures 
were recorded whereas in actuality the primary flow had a 
higher total temperature than the secondary flow. The lower 
temperature reading occurred because an exposed thermocouple 
in supersonic flow measures an adiabatic wall temperature. 

P 

Shapiro introduces the energy equation to describe the 

6 

Shapiro, A. H. , The Dynamics and Thermodynamics of 
Compressible Fluid Flow Volume II (New York: The Ronald 
Press Company, 1954) , pp. 1034-1036. 
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relationship between the adiabatic wall or recovery tempera 
ture the static temperature Tg, and velocity V 

that must exist in the boundary layer over such an exposed 
thermocouple . 




R is the recovery factor that has been shown in many experi 

Q 

ments to be equal to the square root of the Prandtl number 
or 0.85 in laminar flowo Probe II could be used in super- 
sonic flows only if the velocity were known and actual tem- 
perature, the quantity to be determined from measurement. 
This temperature was needed in the ejector tests to deter- 
mine velocity so probe II was not applicable. 

An attempt was made to verify the shape of the pecul- 
iar velocity profiles obtained at axial stations Ml - M3 
by using a hot wire anemometer. The results of such a test 
would be qualitative because of the temperature difference 
between the primary and secondary streams. However, the 
difference in total temperature never exceeded 80°F. With 
the hot wire temperature at 1200 °F, results within 10 per 
cent accuracy could still be obtained. This would be suffi 
ciently accurate to confirm the velocity profile obtained 
with probe I. 

Probe III is shown in Fig. 7. It is a hot wire mount- 
ed on vertical arms attached to 1/4 inch tubing. On each 

9 

Ibid., pp. 1070 - 1072. 
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of six tests the wire maintained its integrity in the secon- 
dary flow but broke as soon as the primary jet was en- 
countered. It is possible that the hot wire broke because 
of vibrations of the thin vertical arms to which the wires 
were attached. 

Probe IV, shown in Fig. 7, was designed such that the 
arms to which the wire was attached were stronger and ori- 
ented parallel to the flow. Both changes were made to 
dampen arm vibration. The probe was placed in the center 
stream before the exhauster was operated. The air flows 
through the turbine and drive nozzle were slowly increased, 
but the wire broke at velocities far below those experi- 
enced diiring normal exhauster operating conditions. 

4 . Description of Test Runs . 

Tests were made to determine the hood total pressure 
variation as a function of drive nozzle total pressure and 
secondary flow rate, to establish the exhauster performance, 
to obtain the variation of mixing pipe static pressure as a 
function of drive nozzle total pressure and secondary flow 
rate, and to study the interaction of the two co-flowing 
streams by obtaining velocity profiles at different axial 
positions of the mixing pipe. Testing was accomplished 
with two different drive nozzles in order to compare the 
effect of variation of the primary flow rate and secondary 
flow annulus area on the problem. 

A total of 34 runs was carried out between February and 
September 1967. Exhauster performance and mixing pipe 
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static pressure rise were obtained for drive nozzle total 
pressure ratios (PTN/^atm^ 2.65, 2.45, 2.18, 1.90, and 
1.63 and for secondary flow rates of 4 and 2 Ibrn/sec. 

The TTTR was set in the following manner. Equation 1 
given in Sec. 5, entitled "Data Reduction", is the equation 
used to determine actual flow rate through the flow measur- 
ing nozzle. Using average values for oc , K, and 



and 


with D 2 


known 


for the TTTR installation, an expression 


for 


the differential 


pressure 




across the pressure taps 


was 


obtained. 
















[."Si 




] 



w = desired secondary flow rate 
p^ = static pressure at upstream tap 

= fluid temperature at upstream tap 

For each test condition, the drive nozzle total pressure was 
set by hand. After allowing time for stabilization of con- 
ditions, readings of the calibrated flow measuring nozzle 
total temperature and pressure p^^ were taken, and an 

h„ calculated that corresponds to the desired flow rate. 

By adjusting the txirbine inlet valve, this value of h^ 
was set. Flow rate regulation to within 0.05 Ibm/sec of 
the desired amount was obtained. 

There was no attempt made to precisely regulate either 

primary or secondary air temperature due to the large number 
of variables and valve settings involved. Turbine testing 
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being carried out concurrent with the ejector tests created 
different turbine efficiencies and operating conditions 
which made temperature variations difficult to regulate. 

The analysis of early test results showed that traverse 
probe position significantly affected hood total pressure 
and secondary flow polytropic efficiencies. The probe was 
removed during final tests to obtain exhauster performance 
and mixing pipe static pressure rise. 

Velocity traverses were initially made at mixing sta- 
tions Ml, M2, M3, M4, m 6, M7 , and M8 of Fig. 3 for a pres- 
sure ratio PTN/P^^ of 2. 18, and secondary flow rates w^ 
of 4, 3, and 2 Ibm/sec. Data points were taken in 1/2 inch 
radial increments at stations M6 - M8. Due to large gradi- 
ents in total temperature and pressure, the increment was 
changed to 1/4 inch for stations M2 - M4 and to 1/10 inch 
at Ml. Additional traverses were obtained at Ml and M4 for 



w 



. of 4 Ibm/sec and PTN/P . of 2.65 and 1.63. 
t atm 



It was not possible to change the traverse probe posi- 
tion without shutting down the exhauster system. After 
restart and stabilization of the system, the new operating 
conditions could vary in total temperatvure up to 3 degrees 
from the previous run. The consequence was considered 
secondary compared to the effect probe position had on ex- 
hauster operation. 

Testing with probe II in run 17 produced total pres- 
sures different from those obtained with probe I for the 
same operating conditions. Investigation of the equipment 
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showed a leak in the traverse probe total pressure tubing. 
The effect of the leak became noticeable with probe II 
because of its smaller total pressure opening. After re- 
pairs, a retest of data points indicated the leak had 
existed during all traverse tests. 

The planned test program was revised and retests ini- 
tiated. Final tests included obtaining the hood total 
pressure and mixing pipe static pressure rise versus drive 
nozzle total pressure, for secondary flow rates of 4 and 2 
Ibm/sec. Velocity traverses were taken at stations Ml - M8 
(less M5) for w^ of 4 and PTN/P^^jjj of 2.18. The tests 
were carried out for both drive nozzles. 

Several attempts to verify the velocity profile at 
station Ml with a hot wire anemometer proved futile. The 
wire could not withstand the impact of the primary stream 
and broke immediately upon insertion in every test. 

5 . Data Reduction . 

Secondary Flow 

Secondary flow rates were measured by a flow nozzle 
between tanks 1 and 2. The basic formula for flow rate in 
Ibm/sec was given by Eckert^*^ to be 



W = 0.1 If 




(1) 



10 



Eckert, op . cit . , Eq. 145, p. 98. 
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D 2 = nozzle throat diameter (in.) 

OC = coefficient of thermal expansion based on 
2 

(dimensionless) 

K = discharge coefficient (dimensionless) 

Y-^ = expansion coefficient (dimensionless) 

p^ = static pressure at upstream tap (in. Hg.) 

= fluid temperature at upstream tap (°R) 

h^ = differential pressure across taps (in. H 2 O) 

Values for OC were determined as outlined by Eckert. 

The expansion factor was changed to that given by the 

ASME Power Test Codes. Discharge coefficient K was 

13 

found by Naviaux to be a function of Reynolds number. 
Naviaux plotted nozzle flow coefficient versus Reynolds 
number and by curve fitting obtained an eighth order poly- 
nomial for K, with Reynolds miitiber the independent varia.- . 
ble. Computation of the polynomial showed that even the 

highest order terms produce large corrections for K. For 
Rg of 11.31 X 10^ the seventh order correction was +1168.857 

^^Ibid. , pp. 98-99. 

12 

Flow Measurement, Supplement to ASME Power Test 
Codes , New York: ASME, 1959, Chap. 4, Part 5, p. 74. 

13 

Naviaux, op. cit., p. 35. 
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and the eighth order term -96.291. Because values of K 
were determined by Naviaux to be in the small range from 
1.00 to 1.05, the polynomial with such large corrections 
was considered inappropriate, and K was varied in steps 
of 0.01 as a function of Rg using Naviaux 's graphical 



results as a reference. 
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The method used for estimating the leak rate through 

the shaft labyrinth was developed by Egli^^ and revised by 
16 

Eckert to apply to the TTTR installation. 
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oc = discharge coefficient for a single 
throttling (dimensionless) 

=0.76 from Egli for present installation 

2 

A = cross sectional area of seal passage (in. ) 



14 



Ibid. , Fig. 19, p. 68. 



15 

Egli, A. "The Leakage of Steam Through Labyrinth 
Seals," Transactions of ASME, V. 57, 1935, pp. 115 - 122. 



16 



Eckert, op. cit. , pp. 103 - 105. 
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Po. = total pressure at seal entrance (psia) 

p = static pressure at seal exit (psia) 

Tu. = total temperature at seal entrance (°R) 

^o 

n = nurnber of throttlings 

The calculation of secondary flow rate was accomplished 
with computer program SECFLO which is shown in Appendix C. 

Velocity 

The velocity profile data reduction was acconplished 

by the computer program VELOCITY. Total temperatxare, total 

pressure and static pressure were meas\ared. The pressure 

ratio was tested for supersonic or subsonic conditions. 

17 

For supersonic flow, the Rayleigh Pitot tube formula was 
used. 




For subsonic conditions, the standard compressible flow 
relation was used. 



y 




17 

Shapiro, A. H. , The Dynamics and Thermodynamics of 
Compressible Fluid Flow Volume I (New York: The Ronald 
Press Company, 1954) , p. 118. 
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In either case, Mach nvunber was solved by iteration and 
converted to velocity with the measured total temperature, 
us ing : 



The computer program used for the reduction is given in 
Appendix C. 

6. Exhauster Performance Prediction . 

An analysis of the exhauster performance to predict 

the hood total pressure, for specified drive nozzle total 

temperature and total pressure, and secondary air mass flow 

rate and total temperature, has been carried out by Vavra^® 

in the TTTR preliminary design calculations, and used by 
19 

Naviaux. However, both analyses assume that the drive 
nozzle is operating at design conditions, and the static 
pressure at the entrance to the mixing section is constant 
in radial direction. Considering the many variables in the 

18 

Vavra, M. H. , Theoretical Evaluation of Exhauster 
System of Transonic Turbine Test Rig (TN65T, U.S. Naval 
Postgraduate School, Monterey, California, 1965) , pp. 16 - 
21 . 

19 

Naviaux, op. cit. , pp. 22 -27. 




( 6 ) 




(7) 
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physical set-up of the TTTR, such as turbine efficiency, 

RPM, secondary air mass flow rate and total temperature, 
and the range of drive nozzle total pressures available, it 
seems unlikely that the nozzle will operate always at de- 
sign conditions. Test data plotted in Fig. 29 indicate a 
static pressure gradient at station Ml, located just 9/16 
inch from the nozzle exit. Therefore, the possibility of 
a discontinuity in static pressure between the nozzle exit 
(area Al) and the secondary flow annulus (area AIT) at the 

start of the mixing section is contained in the present 

20 

development. This hypothesis is supported by Jaumotte 
who indicated that the classic hypothesis that assumes 
uniformity of the static pressure in the initial section of 
mixing is wrong. In the cited reference it is stated also 
that the laws governing the pressures in the two fluids and 
the factors which influence them are not yet well under- 
stood. 

The complete revised analysis of exhauster performance 
is included as Appendix A. Appendix B gives the computer 
program used to calculate prediction maps of exhauster per- 
formance over its operating range. The pertinent steps of 
the analysis are outlined below with the notations corres- 
ponding to those of Fig. 3. 

20 

Engle, M. O. , Some Problems in the Design and 
Operation of Jet Ejectors, Proceedings of the Institution 

of Mechanical Engineers Steam Plant Group (London, England, 
No. 133, 1963), p. 350. 
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Polytropic process relations as well as non-dimensional 

21 

flow functions are used. These are taken from Vavra. For 
an expansion process 




where 

n = 

n = poly tropic exponent 

dT 

= polytropic efficiency ■ 

^ = ratio of specific heats 

w = flow rate ;(lbm/sec) 

T^. = inlet stagnation temperature (°R) 

^o 

2 

A 2 = flow area at station of interest (in. ) 

= inlet stagnation pressure (psia) 

’“o 

P 2 = static pressure at station of interest 
(psia) 

g = universal gravitation constant 
(Ibm-ft/lb-sec^) 

Siibscript 1 refers to inlet conditions 
Subscript 2 refers to station of interest 



21 

Vavra, M. H. , Problems of Fluid Mechanics in Radial 
Turbomachines (Rhode-Saint-Genese, Belgium: Von Karman 

Institute for Fluid Dynamics, 1965) , VKE Co\irse Note 55a, 
pp. 22 - 35. 
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For a compression process 




= static pressure at 1 (psia) 

„ = rVp 

I- ^O-rip) 

For a choked nozzle 



“ChokeJl 




( 10 ) 



Initially, polytropic efficiencies were assvuned to be 
constant over the operating range of the exhauster. The 
determination or estimation of the values used is described 
in Sec. 8 entitled "Determination of Polytropic Exponents." 



Drive Nozzle 



The flow functions from the nozzle entrance to the 
throat ( f n ) / snd from the nozzle entrance to the exit 
( ) . are 




( 11 ) 



( 12 ) 
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The relationship 
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ex ri, 



(13) 



exists , and 



^e.X ^Ti Af 



(14) 



Expressing the exit flow function in terms of a presr- 
sure ratio 
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(15) 



With ^ given by Eq. 14 and PTN an initial condi- 



tion, PIN may be solved by iteration. With PIN, the 

quantities TIN, VIN, and w may then be obtained. 

n 



Secondary Flow 



For the expansion process from the hood to the mixing 
pipe annulus area, there is 




With secondary flow rate and hood total temperature 
known, a value of PIT is assiamed and PTD calculated by 
iteration and matching. For the assumed PIT and correspond- 
ing PTD, the values of TIT and VlT can be calculated. 
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Mixing Pipe 



With the values calculated for the drive nozzle and for 
assumed secondary air flow conditions, the momentum equation 
was used for the region within the mixing pipe to obtain 
conditions at the diffuser entrance. The basic equation is 



Assumptions include uniform velocity profiles, uniform sta- 
tic pressure distribution over each area, and adiabatic 
flow. The term F£ represents the viscous shear force 
imposed by the wall on the it)oving fluid. Estimation of the 
magnitude of this force is outlined by Naviaux. In the 
calculation, the resultant pressure drop produced by 
is determined by a single iteration of Eq. 17 . 

In Eq. 17 , P2 and V2 are unknown. With continuity, 
the pressure P2 can be expressed with the known quanti- 
ties and V2/ which leaves one equation to be solved with 
the unknown quantity V2. With V2/ all conditions at the 
diffuser entrance can be obtained. 




Diffuser 



For the compression in the subsonic diffuser 




( 18 ) 
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The pressxire P 3 is obtained by iteration; and if it does 
not match atmospheric conditions, a new value of PIT for the 
secondary air must be assumed and the calculation repeated. 

7 « Test Results and Discussion . 

Hood Total Pressure 

Tests of the exhauster system were conducted to measure 

the hood total pressure obtained over the operating range of 

the drive nozzle. Runs were made at a secondary flow rate 

of 4 and 2 Ibm/sec for each drive nozzle. 

Two attempts to obtain the information were discarded. 

Rvin 11 was made with the traverse probe at station M2 of 

the mixing pipe. Comparison with other test data showed 

that hood total pressure was increased by as much as 1 inch 

of mercury with the traverse probe installed, particularly 

when the probe was close to the drive nozzle. The increase 

in hood pressure is the result of losses due to shock waves 

that attach to the probe, the drag of the probe in the 

stream, and the wake formed behind the probe. Run 29 with 

drive nozzle II was not used due to the discovery that a 

large seal in the hood had not been properly installed. 

Run 18 and 22 for drive nozzle I and run 32 for drive nozzle 

II were successful. The results are plotted in Fig. 11. 

The graph shows that for the same secondary flow rate, 

drive nozzle I creates a lower hood total pressure than 

nozzle II. This occurs primarily because nozzle I has a 

2 

throat area of 10.42 in. and nozzle II a throat area of 
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8o29 in. 2 so that with choked conditions, about 1.51 Ibm/sec 
greater flow rate of primary air exists with nozzle I. 

Hence, its driving momentum flow is larger. 

For drive nozzle I, a lower hood pressure resulted for 
a smaller secondary flow rate at a given drive nozzle pres- 
sure ratio. In a relative sense, the effect is the same as 
described above. For the same drive nozzle pressxare, which 
establishes the mass flow rate and momentum flow of the 
primary nozzle, it is possible to produce lower pressures 
in the hood for smaller secondary flow rates . 

Mixing Pipe Static Pressure Rise 

Mixing pipe static pressures were recorded for each 
pressure ratio PTN/Pg^^^j^ and flow rate. Graphs showing 
the static pressure variation as a function of distance 
from the drive nozzle exit are given in Figs. 12, 13, and 
14. In general, there occurs a rapid increase in the sta- 
tic pressure along the initial region of the mixing pipe. 

In the main mixing zone where the primary and secondary 
core regions have disappeared, there is a slight but steady 
decrease in static pressure in each case, except for 
PTN/P^^jjj = 2.65 in Fig. 12. The exception could be due to 
the fact that incomplete mixing has occurred in the mixing 
pipe for the extreme condition tested. This point is the 
greatest drive nozzle pressure ratio used and lowest hood 
total pressxare produced during testing. As expected, the 
static pressure rises in the diffuser to the atmospheric 
pressure at the exit. 
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No attempt was made to predict analytically the mixing 
pipe static pressure as a fvinction of axial location. Due 
to the extremely complex nature of the flow which includes 
mixing, turbulence, and shocks in this region, the momentum 
equation used in the exhauster prediction analysis is 
applied only for the whole mixing pipe to obtain average 
conditions ahead of the diffuser. The rapid rise of static 
pressure in the initial region is due to the slowing down 
of the supersonic stream as it mixes with the subsonic 
secondary air. Also present is a series of expansions and 
shock waves to equalize the static pressure between the two 
flows. A possible extension of this thesis would be to 
devise a mixing chamber for Schlieren photographic study 
of the mixing phenomena in compressible flow. 

The slight but steady decrease in static pressure in 
the main mixing zone indicates the loss of available energy 
due to losses caused by the shear forces in the flow and 
the continued turbulent interchange between particles. 

Figures 13 and 14 show that at L/D of about 4, the 
wall static pressure is higher for a drive nozzle total 
pressure ratio of 2.45 than for 2.18. This condition re- 
peated itself on 4 different occasions, but each time the 
manometer board was very xanstable with pressure variations 
at stations Ml M3 up to plus or minus 3 in. Hg. observed 
out of 6 - 8 in. Hg. vacuum. The unsteady condition was 
not reflected in the hood total pressure. It appeared to 
be the result of a shock moving fore and aft in the initial 
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region of the mixing pipe. A very definite sound variation 
was heard while operating at this condition which would 
tend to support the hypothesis of an lonstable shock move- 
ment. 



Velocity Profiles 

Velocity profiles obtained with probe I for drive 
nozzle I are shown in Figs. 15 through 19, and those for 
drive nozzle II in Fig. 20. There occur large gradients 
between the two co-flowing streams at station Ml and con- 
siderable non-uniformities of the velocity profiles. Both 
of these effects decrease with increase in axial distance 
from the drive nozzle. 

The non-uniformity of velocity profile had been ex- 

22 

pec ted because of Naviavix’s results but not to the degree 
observed in the ejector tests. However, Naviaux did not 
traverse closer than 4 inches to the drive nozzle, whereas 
at station Ml the traverse total pressure opening is only 
9/16 inch from the drive nozzle exit plane. 

The xinusual profiles are probably caused by the shock 
structure that actually exists in the flow. A full discus- 
sion with sketches of the possible pattern is included in 
Sec. 10. Essentially, the presence of internal oblique and 
normal shocks causes sharp gradients in pressure and velo<c- 
ity in the flow with accompanied differences in flow direc- 
tions. The large peak in static pressure shown in Fig. 29 

22 

Naviaux, op. cit. , Fig. 3, p. 52. 
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for station Ml is probably due to a large radial velocity 
component in the mixing zone at this radial position. 

Abramovich cites experiments by Zhestkov for a free 
jet operating at off-design conditions. The results are 
similar to the ejector tests with regard to the decrease 
in velocity in the center of the primary jet, compared to 
that in the jet* s boundaries . By comparing the velocity 
profiles at stations Ml, M2, and M3 it is apparent that the 
secondary air increases in velocity and the secondary core 
region decreases in size. The primary stream velocity 
decreases and the primary core area decreases as the mixing 
zone increases in size. 

Interesting results were obtained by non-i 

dimens ionalizing the velocity profiles obtained in the 

main mixing region of the mixing pipe, in the manner given 
24 

by Abramovich for compressible flows. The main mixing 
region is the part of the mixing pipe where the "core" 
flows of primary and secondary air have disappeared. 
Appendix D describes the regions more fully. Abramovich 
uses the ordinate V/V^ where V is the velocity at a given 
location, is the axial velocity. For the absissa. 



I 



y 

2.21 Yc 
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Abramovich, op. cit . , pp. 335 - 338. 
94 

Ibid. , pp. 275 - 276. 
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where y and represent the coordinate of the point 

which corresponds to an arbitrary value of the dimension- 
less velocity, and the coordinate where the velocity is 1/2 
the axial velocity^ respectfully. Figures 21 and 22 graph 
the result. In both cases a single profile is obtained. 

Station M4 is not inside the main mixing region with 
drive nozzle I, but it is with drive nozzle II. Figure 23 
shows that the non-dimensional profile at station M4 for 
nozzle I deviates from the single profile obtained by non- 
dimens ionalizing the velocity profiles of stations M6 - M8. 
From these results, it appears that a sufficient condition 
to determine the start of the main mixing region in a 
restricted jet is to determine where the velocity profile 
may be non-dimens ionalized to a single profile. The 
following paragraphs discuss the mixing regions in greater 
detail. 



Mixing Zone 

The shape of the mixing zone determined from the velo- 
city profiles is shown in Figs. 24 and 25 for nozzle I and 
nozzle II, respectively. This is at best an approximation 
because of the limited number of traverse positions in the 
initial mixing region. The mixing zone appears to touch 
the mixing pipe wall at an L/D of 3.78 for nozzle I, and 
4.71 for nozzle II. D is the nozzle exit diameter. 
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25 

Victorin. carried out a series of experiments using 

water as the driven and driving fluid and determined that 

the velocity distribution in the mixing chamber behaved in 

the same fashion as a circular wake with the width proper - 
1/3 

tional to X , where x is the axial position. The 
assumption is made that the primary fluid emanates at a 
point source at x = 0 located at a position upstream of 
the nozzle so that the velocity distribution would have 
spread to match the nozzle exit dimensions. Calculations 
of the theoretical distance , to predict where the 

mixing zone meets the mixing pipe wall, gave an L/D of 8.27 
for nozzle I, and 9.71 for nozzle II. These conditions are 
shown in Fig. 26. 

It appears from the traverse information that the out- 
ward spread of the mixing zone occurs at a faster rate than 
the incompressible approximation. From the test data ob- 
tained, a closer approximation for the spread of a jet in 

2/5 

compressible flow would be y = x . The relation predicts 
that mixing zone touches the pipe wall at an L/D of 3.87 
for nozzle I, and 4.86 for nozzle II. 

From the traverse data from ;cuns 4-17 it can be con- 
cluded that the shape of the mixing zone is independent of 
secondary air mass flow rate. This conclusion is valid 
even though a leak existed in the traverse probe total 

25victorin, K. , Investigation of Turbulent Mixing 
Processes, NA.CA TM 10 96, 1946. 
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pressiare line dviring these runs. The leak did not affect 
the radial position of the steep total pressure gradient 
observed in test data. The radial position where the steep 
gradient in total pressure occurred was invariant as the 
secondary flow rate was set at 4, 3, and 2 Ibm/sec. 

8. Determination of Poly tropic Exponents. 



The theoretical model which has been assumed for the 
prediction of exhauster performance uses flow functions 
with polytropic exponents to describe the flow in four 
regions of the exhauster; namely, nozzle entrance to nozzle 
throat (n^) , nozzle entrance to nozzle exit (n^^) , secondary 
air from the hood to the annulus at the start of the mixing 
pipe (n^) , and the compression process in the diffuser 
(n^) . The relationship between polytropic efficiencies 



General 



( ) and polytropic exponents (n) , for an expansion pro- 




cess, is 



n, 





and for a compression process 



n 



c 
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Eckert 



assumed polytropic efficiencies to be: 
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7^^ = 0.975, 72, = 0.925, 7^^ = 0.90, and = 0.70. 



Secondary Flow 



Tests were conducted to determine n^j and 72, 
the secondary flow. Secondary flow rate (w^) was maintained 
at 4 Ibm/sec, and drive nozzle total pressure ratio 
PTN/P^^ was set at 2.65, 2.45, 2.18, 1.90, and 1.63. The 
flow fxinction for this expansion is 



Quantities on the left were either measured or are con- 
stants. The pressures on the right were measured, n^j was 
determined by iteration in computer program PDLY. Results 
of the tests are given in Table I. For 4 Ibm/sec secondary 



for the supersonic operating range of the exhauster. Con- 
ditions during the test were stable. 

Initial testing to determine n^ for secondary flow 
rate of 2 Ibm/sec was not immediately successful. In two 
separate tests, results showed polytropic exponents from 
1.08 up to 1.41^ and one value of 1.61 which is impossible 

2 6 

Eckert, op. cit. , p. 93. 




(19) 



flow rate, the average values are n^j = 1.32 and 
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since cannot be larger than 'if =1.41. During the 

test, annulus static pressure was unsteady and an average 
reading was used for data reduction. Figure 13 indicates 
that an unusual condition exists in the mixing pipe between 
pressxire ratios of 2,45 and 2.18, The curves showing mix- 
ing pipe static pressures cross each other for the given 
conditions. This phenomenon repeated itself on four sepa- 
rate occasions. 

In order to determine the effect of a small error in 
the annulus static pressure readings, the computer program 
POLY was used. Values for hood total temperature, total 
pressure, and secondary flow rate were held constant and 
values of annulus static pressure varied in increments of 

0. 1 inch mercury or 0.05 psia. The result is given in 

Table II. It is evident that annulus static pressure is 
critical to the calculation of since a change of 

1. 0.1 psia from the test value changes between 1.39 

and 1.24. In the final test it was attempted to take all 
readings simultaneously. Results of this test are con- 
sidered accurate and are given in Table III. They show 
that n^ is about 1,28, with ^ equal to 0.84. The 
tests showed that the polytropic efficiency was independent 
of drive nozzle total pressure but dependent on secondary 
flow rate. 



Diffuser 

Diffuser polytropic efficiency was determined empiri- 
cally from the static pressure rise between the mixing pipe 
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station M9 to the diffuser exit for all operating conditions 
tested. Figure 27 is a graph of the diffuser static pres- 
sure rise for secondary flow rates of 4 and 2 Ibin/sec for 
different pressure ratios of drive nozzle total pressure and 

atmospheric pressure. 

Initially, the exhauster performance prediction program 
was used with the actual exhauster operating conditions as 
inputs and with a diffuser efficiency 7^^ =0.7 as assumed 
by Eckert. The predicted static pressure rise from M7 to 
atmospheric conditions based on Tl^ of 0.70 is plotted also 

in Fig. 27 for 4 and 2 Ibm/sec secondary flow, showing that 
the diffuser efficiency is larger than 0.7. A diffuser 

efficiency of 0.90 was assumed, and the results of the cal- 
culated static pressure rise in the diffuser is shown in 
Fig. 28. Agreement between calculated and actual diffuser 
static pressure rise is good, but the graph shows that there 
is a small variation of diffuser efficiency with drive noz- 
zle total pressure ratio. A diffuser efficiency of 0.90 
very accurately predicts the performance at the lower drive 
nozzle total pressure ratios up to about 2.1 but under- 
estimates the performance at the higher pressvure ratios. 

The prediction program was altered to vary with drive 

nozzle total pressure. Values used for diffuser efficiency 
at different drive nozzle total pressure ratios are shown 
in Table IV. The results of calculated pressure rise in 
the diffuser compared to test data are in Fig. 28. The 
curves show good agreement between calculated and test 
results . 



4,5 



It was noted also that the diffuser efficiency varied 
slightly with secondary flow rate. This effect is small 
and need not be taken into account in the simplified model 
used for the prediction program. 

9. Determination of Average Shear and Friction Coefficient. 

Mathematical solutions do not exist that exactly pre- 
dict the conditions in a flow as complicated as that in the 
mixing ttibe of the exhauster. However, it was attempted to 
establish a value for the average shear stresses in the 
flow by an approximate method. Such an approximation was 
obtained by not attempting to satisfy the hydrodynamic dif- 
ferential equations for each particle but by integrating 
the profile distributions obtained in laboratory tests and 
hence solving the momentiam equation for the flow between 
two stations. The momentum equation 




where; 1 represents conditions at the upstream station 

2 represents conditions at the downstream station 

Ag = surface area of the mixing pipe between 
stations 

can be used, but in turbulent flow with large velocity gra- 
dients 'T' consists of the "true shear" due to the visco- 
sity of the fluid^ and the "apparent shear" due to the losses 
produced by the turbulent interchange of fluid particles. 



46 



27 

Prandtl states that except in layers directly next to the 
wall, the true shear is small enough to be neglected in com- 
parison to the apparent shear. The approach assxaraes axi- 
symmetric flow, a perfect gas, axial velocities everywhere 
and steady flow, with gravitational effects being ignored. 
Equation 20 is used for the profiles of velocity, pressure 
(total and static) , and temperature measured at stations 
Ml, M4, and M8. 

With the pressures in psi and r in inches, there are 



dA = ZTT rir 

I • au I n IMW -p \/ Ax - 

<km = Py Af\ - V I?, iz ~ g -p / 

Equation 20 becomes 




Z1T 
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Niaraerical integration of large scale reproductions of Figs. 
29 and 30 gave at station Ml 



2TT 





3HZ.y + mr.O 



an.s lbs. 
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Prandtl, L. , The Mechanics of Viscous Fluids, 
Aerodynamic Theory (California: Durand Reprinting 

Committee, California Institute of Technology, January 
1943), Vol. Ill G, p. 135. 
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at M4 



= 33H.H t 5" 2 S'. 5" x 8^3.^ \^s. 

H 



1 -p rJLr - S6S.3 ^ 181.S IW 

Id ^ 

The results, which show an increase in momentum from 
station Ml to M4, must be in error although use of the 
momentum equation to model the mixing pipe performance and 
the mathematical integration can be considered correct. 

The error could have been introduced by the bloclcage effect 
of the traverse probe and the variation of the overall 

exhauster performance with probe position. To illustrate 
this condition. Table V shows the values of the hood total 
pressxire for the same exhauster operating conditions, but 
with the traverse probe at different locations. A varia- 
tion of 1.65 in. Hg. hood total pressure occurs between 
the conditions where the traverse probe is at station Ml, 
and the probe is completely removed. Another possible 
reason for the error is that it was not possible to produce 
exactly the same settings of pressures and temperatures for 
each traverse since, to change traverse probe locations, it 

was necessary to shut down the turbine and exhauster. Re- 
establishment of all previous conditions was extremely 
difficult. Analysis of the data indicates, however, that 
this effect should be secondary. 






( 5 ) 

and at M8 
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The calculation was continued between station M4 and 



M8 using the integrated values obtained. 



Ff = 869. “I - 787.3 



fis 



6.96 



T 



8.19 



82. 6 IW. 



From the test data and results obtained, the turbulent 

friction coefficient f was also calculated. Using 
28 

Vavra, the pressure drop A p due to friction may be 
approximated by 

r L 8 V ^ 

AV = + f-ps-f) 



where 



with 

p = 

■f = 



f = turbulent friction coefficient (dimension- 
less) 

L = length of pipe (in.) 

D = diameter of pipe (in.) 

p = average density between stations (slugs/ft^) 
V = average velocity between stations (ft/sec) 



_ &.38 JP 

RT % lOO T 



179.1 



T A-P 

p v*L 



('sLx/js/ft") 



28 

Vavra, M. H. , Theoretical Evaluation of Exhauster 
System of Transonic Turbine Test Rig (TN65T, U. S. Naval 
Postgraduate School, Monterey, California, 1965) , p. 26. 
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The total frictional force may be expressed in terms of a 
pressure drop 

= ^7^ ^1° =0.303 (lbs.) 

With Ff determined, f becomes 

Average values used to calculate f were; T = 530°R, 
p = 12.4 psia, and V = 700 ft/sec. 

f = 0.106 

This is one order of magnitude larger than the friction 
factor 0.017 predicted by the Moody diagram for commer- 
cial steel pipes. The value should have been closer even 

with velocities of 700 ft/sec. Shapiro cites numerous 
30 

tests that confirm the use of the incompressible friction 
formula for conditions up to Mach = 1.0. 

The calculations were repeated between Ml and M8. 

With Ff = 29.7 lbs., L = 86 in. , p = 11.3 psia, T = 500°R, 
and V = 750 ft/sec, Eq. 21 gives 

f = 0.024 

This value is still large but has the correct order of 
magnitude for the friction coefficient. 

29 

Shapiro, op. cit. , Fig. 6.15, p. 184. 

^°Ibid. , p. 1071. 
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The results of the calculations for friction coeffi- 



cient suggest that the traverse data at station M4 are not 
accurate. Better results were obtained between Ml and M8 
than M4 and M8. No particular quantity in the M4 data was 
considered unreasonable. The conclusion is that the cal- 
culation of an average shear force and friction coefficient 
for the mixing pipe, using the ejector test data, produced 
unreliable results. 

10. Evaluation and Discussion of Predicted Results . 

Actual test conditions were input to the exhauster 
prediction program to determine a predicted hood total 
pressure. Three predictions were made: Program EJECT 1 

used constant polytropic efficiencies determined from test 
data for the secondary air and diffuser and an assumed 
drive nozzle polytropic efficiency of 0.975, EJECT 2 used 
Eckert's polytropic efficiencies, and EJECT 3 used EJECT 1 
efficiencies, except that the diffuser efficiency was 
varied with drive nozzle total pressure from 0.90 to 0.93. 
The results of the predictions for drive nozzle I and a 
comparison with the test results of run 18 are shown in 
Fig. 31 for 4 Ibm/sec secondary flow, and in Fig. 32 for 
2 Ibm/sec. Appropriate changes were made in the program 
to reflect the new dimensions of drive nozzle II. The 
prediction compared to test results is shown ..in Fig. 32. 

The results show that Eckert's assumed polytropic 
efficiencies greatly underestimate the exhauster performance. 
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At high drive nozzle pressure ratios, the use of the theo- 
retical prediction program EJECT I with experimentally- 
determined exponents very closely establishes the actual 
performance of the exhauster. The accuracy of the predic- 
tion decreases as drive nozzle pressure ratio decreases. 

At the low pressure ratios the exhauster performs much 
better than predicted. 

Much of the difference noted above can be attributed 
to the physical formation of an oblique shock within the 
drive nozzle for off-design operation. The classical 
approach to describing nozzle operation notes a third 
critical pressure which is the exit static pressure that 
provides design operation. As the receiver static pres- 
sure is increased, an oblique shock forms, originating 
at the exit plane edge. Further increase in static pressure 
increases the angle of the oblique shock until a normal 
shock forms in the exit plane for the second critical pres- 
sure. Still further increase in pressure causes the normal 

shock to move from the exit toward the nozzle throat. This 

idealized flow development does not occur in actuality 

because of friction and three-dimensional flow effects. 

The latter, in particular, produce: internal oblique shocks 

which have the effect of extending the shock action over an 

extended axial distance in the nozzle. Tests conducted by 
31 

Stodola show internal shock effects over the axial distance 
of a nozzle. 

^^Hall, N. A., Thermodynatmics of Fluid Flow (Englewood 
Cliffs, N. J. : Prentice Hall, Inc., 1951), Fig. 9.5, p. 136. 
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By removing the mixing pipe and operating as a free 
jet, the presence of an internal oblique shock in the drive 
nozzle could be made visible because of condensations of 
the atmospheric humidity in the primary stream. Figure 33 
is a photo of the free jet operation. The oblique shocks 
originated inside the nozzle and extended outside the 
nozzle to their intersection with a near normal circular 
shock in the center area of the nozzle. A sketch of the 
observed pattern is shown in Fig. 34. With increased drive 
nozzle pressure, the entire shock structure moved towards 
the nozzle exit. 

To determine the displacement of this shock structure 
in the drive nozzle under exhauster test conditions, six 
static pressure taps were arranged in drive nozzle I in 
positions shown by Dl through D6 in Fig. 35. At a secon- 
dary flow of 4 Ibm/sec the drive nozzle total pressure was 
set at values of PTN/P^^ of 2.65, 2.45, 2.18, 1.90, and 
1.63 to duplicate conditions of rvin 18. The measured 
static pressures at the nozzle walls are shown in Fig. 36. 
They show that the internal shock moves toward the nozzle 
exit if the pressure ratio is increased. It should be 
noted that for each condition the pressure at the discharge 
of the secondary nozzle was below that required for oblique 
shock formation at the exit plane for isentropic flow 
theory; that is, the pressure was always less than the 2“ 
critical pressure. 
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During the free jet operation it was observed also that 

the primary jet did not completely fill the exit plane when 

oblique shocks occurred inside the nozzle. There existed a 

slowed down region of flow at the outer perimeter of the 

nozzle discharge that had a radial thickness of 0.5 to 0.75 

incho This phenomenon can be noticed in Fig. 34 and must 

be responsible for the apparent thickness of the mixing zone 

at the nozzle exit which is illustrated in Fig. 24. 

The order of magnitude of expected boundary layer 

thickness at the nozzle exit was determined with Eq. 22 

which can be used to evaluate the boundary layer thickness 

32 

on a flat plate in turbulent flow. The Reynolds nurtiber 
applied is based on distance from the nozzle throat to its 
exit. 




^ = boundary layer thickness 
V' = kinematic viscosity (ft/sec 2) 
i = distance from throat to exit (ft) 

V = velocity (ft/sec) 

Based on average conditions 

^ ^ 0.07 inch 

This boundary layer thickness is much smaller than the ob- 
served distance of about 0.5 to 0.75 inch. The increased 

^^Schlichting, H. , Boundary Layer Theory (New York: 
McGraw Hill Book Co. Inc., 1955), Eq. 2.9, p. 38. 
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boundary layer thickness must be caused by flow separations 
from the nozzle walls. This separation is the result of 
the internal oblique shock structure that exists in the 
drive nozzle. The nearer the throat the separation occurs, 
the less of the diverging portion of the nozzle is used to 
increase Mach nvimber of the flow; that is, full e3q>ansion 
is not realized. The effect is readily apparent by ob- 
serving average nozzle Mach numbers obtained at station Ml 
for ^"^/Patm 2.65, 2.18, and 1.63. Test data showed 
Mach numbers of 2.08, 1.18, and 0. 88 respectively. Figure 37 
is a sketch of the shock conditions that appear to exist. 

The effect of the internal oblique shock and resultant 
primary flow separation is to raise the static pressure of 
a portion of the nozzle exit area to the pressure of the 
static pressure at the discharge of the secondary nozzle. 

The momentvun theorem expressed by Eq. 17 uses in the 
actual prediction program 

PSUM A2 - + P^^ A^^ 

A 2 = mixing pipe area 

= secondary flow annulus area 
= nozzle exit area 

= nozzle exit plane static pressure 
= secondary air static pressure. 

In the program the value of PSUM is increased by itera- 
tion until the conditions at the diffuser exit match the 
atmospheric pressure. In the prediction program the 
assumption is made that P^^^j and Pj^ij, are constant, 
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occurs 



respectively, over and A^^. In reality 

on an area smaller than A^, and on an area greater 

than for conditions where the internal oblique shock 

exists. The result is that for a given value of PSUM A2, 
the pressure Pj^^ may be smaller than predicted to satisfy 
the equality. Thus, the hood total pressure would also be 
lower. This condition should be taken care of in the pro- 
gxam for a more correct performance evaluation at low 
pressure ratios. 

11. Conclusions and Recommendations . 

The performance analysis of the exhauster system based 
on polytropic efficiencies and conservation of momentum 
very closely predicts the exhauster performance with re- 
gaurds to hood total pressure at high drive nozzle pressure 
ratios. The prediction under -estimates exhauster perform- 
ance at the low pressure ratios because of the presence of 
internal oblique shocks in the drive nozzle. Accuracy ob- 
tained by the prediction program was the same for drive 
nozzle II which had dimensions different from nozzle I. 

The hood total pressure is significantly effected by 
the presence of a traverse probe in the mixing pipe, ohe 
amount of hood pressure rise is a function of probe posi- 
tion. Therefore, traverse information presented cannot be 
considered as having been obtained for precisely the same 
operating conditions. 
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Traverse profiles near the drive nozzle exhibit steep 
gradients of pressure, temperature, and velocity. Velocity 
profiles obtained appeared irregular. The lower centerline 
velocity is probably due to the shock structure that occxirs 
in the nozzle. It appears that a sufficient condition to 
determine the start of the main mixing region in a restric- 
ted jet is to determine where the velocity profile may be 
non-dimens ionalized to a single profile. Results from flow 
traverses showed that the initial mixing zone elands 
radially at a faster rate than predicted by Vic tor in. s 
incompressible approximation. 

Secondary flow polytropic efficiencies are independent 
of drive nozzle pressxare but dependent on secondary flow 
rate. Diffuser polytropic efficiency is independent of 
flow rate and slightly dependent on drive nozzle pressure. 
The diffuser polytropic efficiency is about 0.90 which is 
considerably higher than the previous estimate of 0.70. 

Little test data are available on the behavior of com- 
pressible free jets. The facility used for the ejector 
tests could be adapted to obtain test data for experiments 
in this area. In addition, further insight into the shock 
pattern in the drive nozzle and the mixing pipe might be 
gained if an adaptation for Schlieren photography were de- 
signed. A- study of the effects of the frictional and 

three-dimensional effects on shock formation in a circular 
nozzle is considered the next step necessary to improve the 
exhauster prediction analysis so that it produces more 
accurate results for the entire operating range. 
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TABLE 



ANNULUS 



POLYTROPIC EXPONENT 



STATIC PRESSURE 
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Actual Test Reading 



VARIATION OF 

Wira ASSUMED ANNULUS STATIC PRESSURE 
Wj = 1.92 Ibm/sec TTD = 547.5® PTD =6.8 PSIA 
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ric 
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Variation of Diffuser Polytropic Efficiency 
With Drive Nozzle Pressure Ratio 



TABLE IV. 
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Hood Total 
Pressure 



Traverse 
Probe Location 
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Hood Total Pressure 
versus 

Traverse Probe Location 
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EJECTOR TESTS 



TABLE V. 
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SHAFT LABYRINTH SEAL CONFIGURATION 
TRANSONIC TURBINE TEST RIG 

FIGURE 4 
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DRIVE NOZZLE I DESIGN 
EJECTOR TEST SERIES 
FIGURE 5 
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TRAVERSE PROBES 
EJECTOR TEST SERIES 
FIGURE 7 
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FIGURE 
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TRAVERSE PROBE INSTALLATION 
EJECTOR TEST SERIES 
FIGURE 10 
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Mixing Pipe 
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VELOCITY PROFILE -STATION Ml 
PTN/f^lm =2.18 Wy =4.0lbm/sec. 

EJECTOR TESTS 
DRIVE NOZZLE I 
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VELOCITY PROFILE - STATION M 4 
PTN/Po|^= 2.18 =4.0 lbw/»®C. 

EJECTOR TESTS 
DRIVE NOZZLE I 
FIGURE 18 
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Mixing Pipe 
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VELOCITY PROFILES 
PTN /Pn4m= 2.48 Wt =4.0 Ibm/sec. 
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MAIN MIXING REGION 
NON-DIMENSIONAL VELOCITY PROFILE 

EXHAUSTER NOZZLE II 
PTN/Pq,^ *2.18 w-r = 4.0lbm/»ec. 
FIGURE 22 
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MIXING STATION M4 
NON-DIMENSIONAL VELOCITY PROFILE 
DRIVE NOZZLE I 
FIGURE 23 
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DRIVE NOZZLE I 
MIXING ZONE 

PTN/fg,tn,=2.l8 Wj-4 Ibm/sec. 
EJECTOR TESTS 
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start of Main Mixing 
Region from Test Data 
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FIGURE 26 
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DIFFUSER STATIC PRESSURE RISE 
COMPARED TO PREDICTION BASE DON EFFICIENCY 7^ =0.70 

FIGURE 27 
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ACTUAL AND CALCULATED 
FIGURE 28 
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EJECTOR TESTS 
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EXHAUSTER PERFORMANCE 

TEST RESULTS AND PREDICTIONS 

DRIVE NOZZLE I 
Wt = 4lbmAec. 
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EXHAUSTER PERFORMANCE 
TEST RESULTS AND PREDICTIONS 
DRIVE NOZZLE IT 
Wy = 4 Ibm/sec. 

FIGURE 33 
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FREE JET OPERATION 
DRIVE NOZZLE I 
FIGURE 34 
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DETAIL OF PRESSURE TAP 
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Station 
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Shock Structure 



APPENDIX A 



Exhauster Performance Analysis 

An analysis was made of the exhauster performance to 
predict hood total pressure which essentially parallels the 
preliminary design calculations of Vavra. However, the 
previous work assumed that the nozzle operates at design 
conditions and that the static pressure at the entrance to 
the mixing section is constant with radius. The analysis 
explained in the following permits a discontinuity in sta- 
tic pressure between the nozzle exit (Area Al) and the 
secondary flow annulus (Area AIT) at the start of the mix- 
ing pipe. 

In order to simplify the equations, polytropic pro- 
cess relations as well as non-dimensional flow functions 
were used. These were taken from Vavra. 

For an expansion process 



Pn / 3 




(Al) 



where 



n = 



y - yo fv-/) 

Lp 



(A2) 
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dT 

dT. 

IS 



n = poly tropic exponent 

= polytropic efficiency | 
^ = ratio of specific heats 
w = flow rate (Ibrn/sec) 



T*. = inlet stagnation temperature (°R) 

o 



2 , 



A 2 = Flow area at station of interest (in. ) 



P. = Inlet stagnation pressure (psia) 

’^o 

P 2 = static pressure at station of interest 
(psia) 



g = universal gravitational constant 
Subscript 1 refers to inlet conditions 



f Ibm ft 
I ^^f sec ^ ' 



Subscript 2 refers to station of interest 
For a compression process 




'y-i 



r (f r- (f r 






(A3) 






static pressure at 1 



(A4) 



The flow through the nozzle will be choked when it is 
operating at supersonic conditions. Therefore $,3 (or 

) must be a maximum value. Differentiating and 
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setting the result equal to zero yield the critical pres- 
sure ratio for choked conditions 






-] 

4rit/c<xl 




Y^- I 



Substituting this value into Eq. Al^ the expression for the 
flow function of a choked nozzle simplifies to 

I 






cWo 

Nonsle 




(A5) 



It was assumed that the drive nozzle total temperature 
and total pressure were known^as well as the secondary flow 
rate and total temperature. Polytropic efficiencies for 
the drive nozzle are assumed or determined from testing and 
held constant over the operating range of the exhauster. 
Notation in the development is consistent with Fig. 3. 



Drive Nozzle 



For a choked nozzle. 




(A6) 




The flow functions from the nozzle entrance to the throat 
( $„ ), and from the nozzle entrance to the exit ( ) are, 

respectively. 
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n 




j 



ex 




There are 

= LA, 



and 




(A7) 



(A8) 



(A9) 



(AlO) 



Expressing the exit flow function in terms of the pressure 
ratio Piju/Piij 




With given by Eq. AlO, and with PTN known, PIN can be 



solved by 


iteration. Then 




Tn 


- fe ) 


(A12) 




= V2iJcpft„-T„y 


(A13) 


w 

^ N 


= 5 P A 


(A14) 
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Secondary Flow 



Secondary air is that part of the high pressure air 
provided by the Allis -Chalmers compressor that passes 
through the test turbine to the hood. For the expansion 
process from the hood to the annulus area 




(A15) 



With the secondary flow rate and hood total temperatures 
known, a value of PIT was assumed, and PTD is calculated 
by iteration and matching. For the assumed PIT 



T. = T, 



TO 



P,r 



ro 






(A16) 



V,^ = V 29 J=p^T,„-T„) 



(A17) 



Mixing Pipe 

The momentum equation was used for the mixing pipe to 
obtain conditions at the diffuser entrance. The flow rate 
through the drive nozzle is known from preceding calcula- 
tions and values must be assumed for the secondary flow 
conditions. Applying the momentum equation to the control 
volume between the nozzle exit and annulus plane (sta- 
tion 1)^ and the entrance to the diffuser (station 2)^ 
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(A18) 




V - ^V - sPA, 

V ''it ^it It 



+ P A - 

'in I 




Equation A18 holds for unlforin velocity profiles and uniform 
static pressure distribution over each area. 

Ihe friction force Ff may be e:q>ressed as a pressure 

loss 

'■f “ T 

where 








L = length of duct 
D 2 = duct diameter 

V = average velocity 
p = average density 
f = friction factor 



(A20) 



The friction factor was determined to be 0.017^ in the man- 

33 

ner outlined by Naviaux. The friction force was figured 
with one iteration only. For the initial solution of 
Eq. A18, Ff was considered zero. After obtaining the con- 
ditions at the diffuser entrance, average values of V and 

p for the mixing pipe were calculated , and a A P found 
using Eq. A20. /^P A 2 was then subtracted from + 

PljjAiN to obtain a new solution of Eq. A18. 

^^Naviaux, op. clt. , pp. 25 - 26. 
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Unknown quantities in Eq. A18 are V 2 and P 2 » The 
pressure P 2 can be e3q>ressed by V 2 to obtain one inde- 
pendent variable^ if the total temperature after mixing is 
taken as the mass weighted average of 







(A21) 



This condition assvimes the existence of an adiabatic pro- 
cess. With 



w„ + w„ = w 
N T 



(A22) 



and the continuity equation 



W - 



s 



(A23) 



V/ _ w 



_ w R 



T - ^ 



(A24) 



Therefore 



■p _ w R 

'2. ~ 



L 



T. 



T^ 



-ji* ] 
^ 3 ^ Cp J 



(A25) 



Substituting Eq. A25 into Eq. A18 

~ iii\/ - P /i +P A 



_ wR I 

' V;, 1. 



■ 






(A26) 
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Multiply Eq. A26 by V 2 and rearrange 

-^WRTTi = 0 



- 

jv wR 


z 

w 


3 





(A27) 



Let 



Fl'- 


ZJcp J 


(A28) 


-h 

J 

1 

II 

CM 

CQ 


^IT * 


(A29) 


B 3 — \A/ 




(A30) 


®1 ’'2 + »2 V2 + B3 = 


0 


(A31) 



Thus 



Ye = 



-B. -'VB/-Me,e. 



I - - _ 






2 6 . 



2 6, 



£ 

.26, 



^\ _ £1 
i. / B. 



(A3 2) 



and 



«?. - n, I if* 



(A33) 



(A34) 



The negative sign before the square root yields the 
stibsonic solution. Supersonic conditions at the diffuser 
are outside to operating range of the TTTR installation. 
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Finally, for the compression in the siabsonic diffuser 



is solved by iteration and compared to atmospheric 
pressure. If it does not match for the chosen values of 
secondary flow conditions, a new PIT is selected and the 
process repeated vintil p^ equals Patm* With P 3 
obtained 




(A3 5) 







(A3 6 ) 




(A37) 




(A38) 
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APPENDIX B 



Exhauster Performance Prediction Program 

The analysis of exhauster performance presented in 
detail in Appendix A was programed in Fortran IV for cal- 
culation on the IBM 360 computer. Input data include the 
assumed polytropic efficiencies, drive nozzle total pres- 
sure, drive nozzle total temperature, secondary air total 
temperature, secondary flow rate and atmospheric pressure. 
The principal quantity to be determined is the predicted 
hood total pressure. The printed output does include 
final values of temperature, velocity and pressure calcu- 
lated during the computations. 

The program uses an executive MAIN routine to arrange 
the sequence of calculations^ but all calculations are 
performed in subroutines. The name and function of each 
subroutine is: 

CSTANT; a grouping of all the basic quantities 
used in the entire program. 

CHOPMI; a calculation of the non-dimensional flow 
fvinction for a choked nozzle. 

NOZZLE; computes conditions at the nozzle exit. 

SECFLO; computes conditions at the secondary flow 
annulus area at the entrance to the mixing pipe. 

MMENIM? Using conditions determined by NOZZLE and 
SECFLO, conditions at the diffuser entrance are calculated 
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based on conservation of momentum. One iteration is in- 



cluded to determine the wall frictional shear force, 

SBCOMP; computes the compression process in the 
diffuser to obtain conditions at the diffuser exit. 

Input information required to operate the prediction 



program is: 
Card 1; 




ETAD - 


Secondary flow polytropic efficiency 


ETAN - 


■ Nozzle polytropic efficiency, entrance to 
throat 


ETAl - 


Nozzle polytropic efficiency, entrance to 
exit 


NSETS - 


' N\jmber of sets of data to be computed 


Card 2 ; 

PTN 


Drive nozzle total pressure (psia) 


TTN 


Drive nozzle total temperature (°R) 


TTD 


■ Hood (secondary air) total temperature (°R) 


WT 


Secondary flow rate 


PA 


■ Atmospheric pressure 



Program EJECTOR with a sample output sheet is given on the 
following pages. 
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SUBROUTINE CHOPHl ( PHI N» ENN, C2 ) 

CHOPHI COMPUTES THE FLOW FUNCTION FOR A CHOCKED NOZZLE 

PHIN=(2./{FNN+l. I )=i'*(l./{ENN-l .) ) *SOR T ( C 2« I ENN- 1 . ) / ( ENN+ 1 . ) ) 

RETURN 

END 



f:XPl=?./ENl 
tXP2={ ENl + 1. I/EWl 
EXP3 = ( ENl- 1. )/ENl 
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APPENDIX C 



Data Reduction Programs 

Data obtained from laboratory tests to compute secon- 
dary flow rate, velocity profiles, and polytropic es^onents 
were reduced by the computer programs included in this 
section. Input requirements are listed on comment cards at 
the start of each program. 

Program SECFLO computes the mass of air that passes 
from the hood to the mixing pipe. The program basically 
calculates Eq. 1 of the main text. In addition, the small 
quantity of air that enters the hood through the shaft 
labyrinth is approximated by Eq. 2. Naviaux's flow nozzle 
coefficient is determined as a function of Reynolds number 
by a series of "IF" statements. The values used range from 
1.00 - 1.05 in increments of .01. 

Program VELOCITY computes Mach number and Velocity 
from traverse probe data. Input pressxires are tested for 
supersonic or subsonic conditions. Either Eq. 4 or Eq. 5 
is used depending on the condition that exists. 

Program POLY computes by iteration the secondary flow 
poly tropic exponent from test data. Eq. 19 with all quan- 
tities known but n^ was used in the calculation. 
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Program VELOCITY Continued 
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APPENDIX D 



Restricted Jet Operation 

In order that the reader might better understand the 
general terms used in this thesis to describe jet operation 
and associated mixing phenomena, the following general 
description of jet operation is presented. 

Figure Dl of this appendix pictures the physical sit- 
uation to be discussed. Essentially, two co-flowing 

33 

streams enter a mixing pipe. Abramovich states that with 
this situation, the velocity components along the y-axis 
(the radial components for the exhauster) prove so small 
that they can be disregarded in engineering problems in- 
volving jet theory. 

The two co-flowing streams enter the mixing pipe with 
their fluid particles traveling at different velocities. 

The viscous forces acting between the two flows cause a 
slowing down of the particles of the primary jet and a 
speeding up of the secondary flow particles. The area of 
the flow where the interaction between the particles of 
the two streams takes place is called the mixing zone. 

The mixing zone thickens with axial distance, causing a 
gradual "eating up" of the primary and secondary core 
regions. 

Two main mixing regions are designated in the mixing 
pipe. The initial region has the axial distance where the 

33 

Abramovich, op. cit. , p. 5 
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secondary core region and the primary core region still 
exist. The main mixing region is that part of the pipe 
where the mixing zone has touched the pipe walls and the 
primary core region no longer exists. The regions are 
shown in Fig. Dl. In the literature use is sometimes made 
of a transition zone which is the region that occurs if the 
mixing zone contacts the pipe wall before the primary core 
is "eaten up" or vice versa. This thesis uses the simpli- 
fied model in which it is assumed that the length of the 
transitional region is equal to zero. In the main mixing 
region, the jet becomes similar in appearance to a flow 
of fluid from a source of infinitely small thickness (in 
an axially symmetric case, the source is a point). 
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